Draft version February 2, 2008 

Preprint typeset using 1^1^^ style cmulatcapj v. ^jllj^A, 



in 
o 
o 

00 



> 

o 
in 
o 

:^ 

Oh 
6 



X 



STATISTICS OF QUASARS MULTIPLY IMAGED BY GALAXY CLUSTERS 
Joseph F. Hennawi^-^-', Neal Dalal^"' Paul Bode^ 

Draft version February 2, 2008 

ABSTRACT 

We compute the expected number of quasars multiply imaged by cluster size dark halos for current 
wide field quasar surveys by carrying out a large ensemble of ray tracing simulations through clusters 
from a cosmological N-body simulation of the ACDM cosmology. Our calculation predicts ^ 4 quasar 
lenses with splittings > 10" in the SDSS spectroscopic quasar sample, consistent with the recent 
discovery of the wide separation lens SDSSJ1004+4112 which has 6 = 14". 6. The SDSS faint photo- 
metric quasar survey will contain ^ 12 multiply imaged quasars with splittings 6 > 10". Of these, 
~ 2 will be lenses with separations 9 > 30", and ^ 2 will be at high rcdshift (zg ^ 4). 

Subject headings: 



1. INTRODUCTION 

Strong gravitational lensing by clusters of galaxies pro- 
vides a unique laboratory for studying the dark matter 
distribution in the largest collapsed structures in the uni- 
verse. The most common manifestation of cluster strong 
lensing are giant arcs, which are highly elongated and 
sometimes multiply imaged background galaxies which 
have been distorted by the deep gravitational potential 
of a massive galaxy cluster. Background quasars can also 
be multiply imaged by foreground lenses. However, un- 
til recently , all o f the roughly 70 lensed quasars known 
ijKochanekl I2004j) had small < 10" separations, ex- 
plained in terms of galaxy-scale concentrations of bary- 
onic matter, despite a number of searches for lenses 
with wider splitting s llMao z et aI][T997: Ofek et al. 2001; 
Phillips et al .1200 It Zhdano v fc Siirdeil2001blMiller et all 
2004HMarble et alJi2004) . 



The discovery of the quadruply imaged quasar 
SDSSJ1GG4+4112 ffnada et a,].l l2nM: 'Og uri et a,lJl2004t 
Unada et al. 2005) in the Sloan Digital Sky Survey 
(SDSS), with a largest image separation of 6 — 14". 6, 
constitutes the first quasar multiply imaged by a clus- 
ter where dark rn atter dominates the lensing potential 
ijOguri et al.ll200^ . Quasars lensed by clusters are much 
rarer than giant arcs simply because the number den- 
sity of quasars behind clusters is much lower than that 
of galaxies. Because massive clusters capable of strong 
lensing are also very rare, large samples of quasars are 
required to find quasars multiply imaged by clusters. 

The advent of large area di gital im aging and spectro- 
scopic surveys, like the SDSS ijYork e t al. 200 3) and th e 
Two Degree Field Quasar Survev l|Croom et alJl200^ . 
has dramatically increased the number of known quasars. 
Although a spectroscopic survey of ~ 30,000 quasars was 
required to disc over one multiply imaged quasar lensed 
by a cluster (fhiada et alJl200,1 lOguri et alJl2004l ). new 
photometric quasar selection techniques applied to the 
current SDSS imaging data promise to yield as many as 
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a million quasars ('Rich ards et al .^'2004V In addition, fu- 
ture large area synoptic surveys such as Pan-STARRS 
and the LSST will be able to use variability selection 
techniques to detect as many as one hundred million 
quasars. The number of quasars lensed by clusters dis- 
covered in these current and future samples of quasars 
will allow statistical studies of the abundance of these 
large separation quasar lenses. 

It has been suggested that the abundance of clus- 
ter lenses can be us ed as a cosmological probe. 
iBartelmann et aP (|1998) found that the predicted num- 
ber of giant arcs varies by order s of magnitude among dif- 
fere nt cosmological models, a nd lBartelmann et al.l l|2003f) 
and iMeneghetti et alJ 1)20041) recently proposed that gi- 
ant arc statistics could be used to distinguish between 
different dark energy cosmological models. Despite the 
fact that the abundance of quasars lensed by clusters will 
be significantly lower, because quasars are point sources, 
the selection function of multiply imaged quasars is much 
easier to quantify than that of giant arcs, because the 
ability to detect giant arcs is extremely sensitive to the 
seeing. Indeed, beginning with Kochanek (1996), the 
abundance of multiply imaged quasars has been recog- 
nized as a potentially powerful cosmological probe. How- 
ever, the dominant uncertainty has always been the accu- 
racy with which the galaxy population can be modeled. 
Because on cluster scales only dark matter and gravity 
are involved, quasars multiply imaged by clusters would 
not suffer from these modeling uncertainties. 

The measurement of a time delay in multiply imaged 
quasar cluster lens could prove to be promising method 
for measuring the Hubble constant, if the cluster lens also 
shows giant arcs. This is because there would be more 
constraints on the mass model than are present in galaxy 
lenses, and it precisely this modeling uncerta inty which 
limits the precision of this technique (see e.g. lKochanei3 
|2003: Kochanek & Schechter 200|). 

Finally, because multiply imaged quasar lenses are 
found by searching for pairs of quasars at the same 
redshift (that is, by surveying the source population), 
whereas giant arcs are discovered by looking for elon- 
gated structures in the vicinity of rich galaxy clusters 
(that is, by surveying the lens population), the two types 
of cluster lens surveys are complementary. 

Nearly all previous investigations of the abun- 
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dance of quasars lensed by massive dark halos 
have used simple spherically symmetric analyti- 
cal profiles for the mass distribution of clusters to 
predict the abundance o f wide separation q uasar 
lenses fMaoz et all 119971: 'Sarbu. Rusin. & Ma 2001!; 
Keeton fc Madau 120011: Wvithe. Tu r ner^fc_Sp_er^ 
2001 : iLi fc Ostrikej 12002 lRusinn200a~lLi fc Ostrike: 
20031 IQguril 120031 IHuterer fc Mai 120041 ILopes fc Mille: 



2004 : IKuhlen. Keeton. fc Madauf 120041 lOeuri et all 
20041 iChenI \200i). However, in the context of gi- 
ant arc statistics, full ray tracing computations 



through clusters from N-body simulations have un- 
equivocally demonstrated that analytical models fail 
to accurately represent strong lensing by clusters 
(jBartelmann. Steinmetz. fc WeissI 11995(1 , underestimat- 
ing the lensing cross sect ions b y as much as two orders 
of magnitude l| Meneghctti. Ba rtelmann. fc Moscardinil 
2003al iHennawi et alJ i2005a,) . In recent progress, 
Qgurifcl^eetonr l|2004|) computed the abundance of 



quasars multiply imaged by clusters analytically, using 
triaxial rather than spherical halos; they found that 
triaxiality significantly increases the expected number of 
lenses. In addition, all previous predictions for the num- 
ber of wide separation multiply imaged quasar lenses 
have neglected the effect of the brightest cluster galaxies 
on the strong lensing cross section, focusing only on the 
effects of the dark matter. However, brightest cluster 
galaxies can en hance strong lensing cross se ctions by 
up to ^ 50% llMmegh^liti^^artelmann. fc Moscardini 
ppOSF, 'Dalai. Holder, fc H ennawil 12004^. 'Ho fc White 
12(30 4; Hennawi et al. 2005a) for small image separations 
{0 < 20"). 

In this paper we compute the expected number of 
quasars multiply imaged by cluster size dark halos in 
current wide field quasar surveys by carrying out a large 
ensemble of ray tracing simulations through clusters from 
a large cosmological N-body simulation. In §13 we briefly 
summarize the ray tracing simulations and our technique 
for adding BCGs to dark matter halos. The details of the 
quasar luminosity functions for the surveys we consider 
are discussed in ^ We present results on the abundance 
of quasars lensed by clusters in § 0] and we conclude in 

2. RAY TRACING SIMULATIONS 

In this section, we briefly summarize the essential 
elements of our strong le nsing simulation s . W e re- 
fer the reader t o [Palal. Holder, fc Hennawil l(2004f l and 
IHennawi et al.l l|2005a|) for details. We will discuss the 
N-body simulations, ray tracing simulations, and our 
method for adding brightest cluster galaxies in turn. 

In order to simulate the strong lensing effects of 
a Universe filled with dark matter, we used clusters 
from a cosmological N-body simulatio n, performed with 
the IV ee-Particle-Mesh (TPM) code of lBode fc Ostrikeil 
(|2003|) . The following cosmological parameters were 
used for this simulation: fl^ = 0.3, J^a = 0.7,/i = 
70 km s~^ Mpc^^, as— 0.95, and ns = l, which are con- 
sistent (within la) of the WMAP derived cosmologi- 
cal parameters (|Spergel et al.l 12003(1 . The simulations 
were performed in a box with a comoving side length 
of L = 320 h-^ Mpc and N = 1024^ particles, giving 
a particle mass of rup = 2.54 x 10^h~^ Mq. The cu- 
bic spline softening length was set to e = 3.2h^^ kpc. 



We use outputs at seven different 'snapshots', covering 
the range of redshifts over which the critical density is 
low enough to produce an appreciable amount of strong 
lensing: z = 0.17,0.29,0.41,0.55,0.70,0.87,1.05,1.26, 
and 1.49. This same simulation ha s been used to in- 
vestigate the statistics of gia nt arcs (|Wambsganss et all 
l2004al IHennawi et al.l 2005a|l . strong lensing cosmogra- 
phv l|Dalal et al .1 1200511 . and the lensing effects of sec- 
ondary matter in the light cone (Wambsganss et all 
I2004htl . 

A "friends-of- friends" (FOF) group finder 
l|Davis et al.l 1198511 with the canonical linking length 
of & = 0.2 was applied to each particle distribution to 
identify cluster size dark matter halos. For each cluster 
with a FOF group mass above MpoF > lO^'^h'^ M©, 
all the particles within a 5h~^ Mpc sphere about the 
center of mass were dumped to separate files and used 
as inputs to our ray tracing code. 

To compute reliable mean cross sections for a clus- 
ter, we must average over many projections to ap- 
propriatel y sample the di s tribut io n of cross s e ctions 
toalal, H older fc Hennawil l200l IHo fc Whitel l200l 
IHennawi et al.ll2005a|l . WfT computed mean lensing cross 
sections by averaging over 125 orientations for all clus- 
ters with Mfof > lO^''-" h-^ Mq, 31 orientations for 
clusters in the mass interval lO^"*'^ Mq < MpoF < 
2Q15.0 j^-i -|y[^^ 3 orientations for all clusters in the 
range lO^* h'^ Mg < MpoF < lO"'^ h'^ Mq. 

The surface density of each cluster projection is com- 
puted on a grid and a deflection angle map is computed 
using FFT methods. Using linked lists, a mapping is 
constructed between source plane pixels and image plane 
pixels. We then compute the statistics of multiply im- 
aged quasars using a Monte Carlo approach. Single pixel 
sources are randomly placed in the source plane. If there 
is more than one image plane pixel corresponding to this 
source plane pixel we count this as a multiply imaged 
quasar. If the image multiplicity is larger than two, the 
separation of the images is taken to be the largest image 
separation in the system, and this is considered a sin- 
gle lensing event, thus we properly count triples, quads, 
etc. The magnifications of each image are also computed, 
which is just the ratio of the number of pixels in each im- 
age to the original number in the source (which is unity). 

Given this list of image splittings and magnifications, 
the number of multiply imaged quasars behind each 
cluster orientation can be computed, given the number 
counts of background quasars as a function of apparent 
magnitude. Our calculation completely accounts for the 
magnification bias, because we use the magnification of 
each image to determine the corresponding number den- 
sity of background sources for the particular survey under 
consideration. 

The ray trace and Monte Carlo must be repeated for 
each projection through the cluster and each source plane 
considered. We ray trace three source planes at Zg = 
1.0, 2.0, and 4.0. The number counts of quasars are 
thus collapsed into three source redshift bins centered on 
these source planes (see given by [0.75, 1.5], [1.5, 3.0], 
and [3.0,5.0], respectively. Because the critical density 
for strong lensing is a slowly varying function of source 
redshift, this binning will not introduce significant errors 
in our calculation. 

Brightest cluster galaxies (BCGs) can en- 
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Fig. 1. — Left: Redshift distribution of quasars for a flux limit of i < 22.7. The green, red, and blue shaded regions indicate the redshift 
bins used to determine the number of quasars for each of the three source planes at Zb = 1, 2, and 4, respectively. The area shaded gray 
represents the number of quasars excluded by our assumption of 60% completeness over the redshift range 2.3 < z < 3.6, where color 
selected samples are incomplete because quasar colors cross the stellar locus. Right: Cumulative number magnitude counts of quasars 
used for each of our three source planes. These curves were obtained by integrating the redshift distribution of quasars over redshift bins 
illustrated by the shading in the left panel. 



hance strong lensing cro ss^ s ections by 



' 50% 
l2003b : 



([Mencghctti. Bartclmann, & Moscardini . 

[Dalai. Holder, fc Hcnnawi .2004 .Ho fc White! 12004 ; 
IHennawi et al.1 2005a) for multiply imaged quasar sepa- 
rations < 20", small enough such that the mass enclosed 
by the critical lines has a significant baryonic com p onent 
Jbalal. Holder. & Hennawi 200^ IHo & Whitel l2004t 
iHennawi et al... ,2005a) . Indeed, the recently discov- 
ered wide separ ation quadruply i maged quasar lens 
SDSS1004-H4112 ijlnada et al.ll2003|) is centered on the 
brightest galaxy of the cluster. Models of the lens which 
include the central galaxy embedded in a dark matter 
halo give an Einstein radius of ~ 7", with the galaxy 
component contributi ng a signific ant amount to the 
lensing potential ( Qguri et al.ll2004D . 

We artificially add baryons, with a singular isother- 
mal sphere profile, to the centers of each cluster by 
'painting' BCGs onto the dark matter surface den- 
sity. Because varying the total mass of the central 
galaxy has a negligible effect on the lensing cross section 
|Dalal. Holder . & Hennawi 2004), we use a fixed frac- 
tion, Mbaryon = 0.003 A/F0F7 of the mass of the dark 
halo. We use a simple scaling prescription for assigning 
velocity dispersions to the BCGs 



300 km s" 



FOF 



2/15 



(1) 



which is based on an empirical scaling relation between 
the optical luminosity o f BCGs and the X-ray tempera- 
tures of galaxy clusters [Edge fc Stewart! l|1991f) and the 
fact that BCGs lie o n the fundamental plane of ell iptical 

;alaxie s ijOeeerle fc Hocssel 199ll) . See .Hennawi et al.! 

20053(1 for details. 

3. MODELING THE QUASAR LUMINOSITY FUNCTION 

We compute the expected abundance of wide sep- 
aration multiply imaged quas ars for two surveys : 
the SDSS spectroscopic survey ((Schneider et al.! !2005(1 
and the fainter SPS S photometric quasar sample of 
!Richards et"an pOOl . 



At low redshift z < 2.3, the quasar lumi nositv func- 
tion h as been measured by several groups fBovl e et alJ 
l200nt !Croom et al.l I200I !Richards et all .2005il . We 
use the dou ble po wer law B-band luminosity function 
l|Bovle et al.!l2000l) 



$(Mb,z) 



lQ0Aif3i + l)[MB-M^iz)] _^ iQ0A(f3H + l)[MB-M^iz)] ' 

where /3/ = -1.64, Ph = -3.43, and = 360(/i/0.50)3 
Gpc^'^mag^i. The evolution of the break luminosity 
M*{z) follows 

M* =M*(0)-2.5(fciz + fc22'), (3) 

with fci = 1.36, fc2 = -0.27,and M^(0) = -21.15 + 
5 log h. 

The quasar luminosity functio n is poorly constrained 
between redshifts 2.3 < z < 3.6. !Wvithe fc Loebl l|2002l) 
used a simple analytical model to fit the do uble power 
law lu minosity function in eqn. Q to both the 'Fan et alJ 
(!2001[1 high redshift ( z > 3.6) luminosity function and 
the !Bovle et al.! l|2000|) low redshift (z < 2.3) luminosity 
function, using a mode l for the evolution of th e break 
luminosity proposed bv !Madau. Haardt. fc Ree3 (p399) 



L^z) = LU0){1 + z 



(„_i)e^^(l + e«^ 



(4) 



where a = —0.5 is the slope of a power law spectral 
energy distribution assumed for the quasar spectrum 
f,j ^ i^" . At the bright end they used the !Fan et all 
l)200l!l slope /3h = 2.58 and assumed a faint end slope 
/3i = 1.64. The other parameters of their fit are = 624 
Gpc-^mag-i, M|(0) = -22.46, = 1.60, ^ = 2.65, 
and ^ = 3.30. This prescription for the luminosity func- 
tion has also been emp loyed in strong lensing studies by 
!Oguri et al.! l|200^ and !Oguri fc Keeto n' f2004*). For red- 
shifts z < 2.3 we use the !Bovle et"aLr |2000) expression m 
eqn. |j2Jl. In the range 2.3 < z < 3.6 we lin early interpo- 
late between eqn. M and the the l!Wvithe fc Loeb 200^ 
fit, and we use the ifWvithe fc Loeb!!2002() fit for z > 3.6. 
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Fig. 2. — Cumulative distribution of image splittings for two diff erent quasar surveys. The left panel is for the SDSS spectroscopic 
quasars and the right panel is for SDSS photometric quasar sample of'R ichards et all 120041) . The green (dotted), blue (dot-dashed), and 
red (dashed) curves are the separate contributions from the source planes at Zb = 1.0, 2.0, and 4.0, respectively. The sum of these three 
curves gives the total number of multiply imaged quasars which is the black (solid) curve. The thin solid (black) line shows the total 
number of lenses if BCGs are neglected and we ray trace through dark matter only. 
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Fig. 3. — Redshift distribution of cluster lenses with splittings 9 > 10" for the spectroscopic quasars (left) and photometric quasars 
(right). The green (dotted), blue (dot-dashed), and red (dashed) curves are the separate contributions from the source planes at 2s = 1-0, 
2.0, & 4.0, respectively. The sum of these three curves gives the total number of multiply imaged quasars which is the black (solid) curve. 



The luminosity functions we are considering are ex- 
pressed in terms of rest frame B-band absolute magni- 
tude, Mb, whereas the SDSS spectroscopic and photo- 
metric surveys have flux limits in the i-band. We thus 
require the cross filter K-correction KBi{z), between ap- 
parent magnitud e i and absolute magnitude B (see e.g. 
IHngg et a].ll2nf)l 

M^^i~DM{z)^KB^{z) (5) 

where DM{z) is the distance modulus. We compute 
KBi(z) from the SPSS c omposite quasar spectrum of 
iVanden Berk et all l|2001f) and the Johnson-B and SDSS 
i filter curves. 

Between redshifts 2.3 < z < 3.6, quasar surveys 
which rely on color selection suffer from incompleteness 
because quasar colors cross the stellar locus (see e.g. 
iRichards etHII I200H IVanden Berk et all l2005j) in this 
redshift range. Because the SDSS spectroscopic and 
photometric surveys rely on color selection, we conser- 
vatively assume that only 60% of the quasars in this in- 
terval are recovered, accounting for this incompleteness. 



An area of ~ 5000 deg^ is assumed for the SDSS spec- 
troscopic quasar survey, which has flux limits i < 19.1, 
for low redshift quasars (z < 2.5) and i < 20.2 fo r high 
redshift quasars {z > 2.5) ijSchneider et al.ll2005(l . We 
assume the faint photometric quasar survey covers the 
larger SDSS imaging area of ^ 8000 deg^, with flux 
limits i < 21.0 for z < 2.5, and i < 20.5 for z > 2.5 
(Richards et al.ll200# . 

Searching for quasar lenses requires requires follow-up 
observations. This is the case for the SDSS spectroscopic 
sample because of fiber collisions^, and also for the faint 
photometric sample because spectroscopy is required to 
confirm that two quasars in a pair are quasars at the 
same redshift, rather than a projected pair of quasars 
or stellar contaminants. The flux limit of the follow-up 
observations is set by the follow-up instrument and is 
typically fainter than the flux limit of the parent sample. 
Thus for both surveys, we assume a flux limit of i < 21 

^ The finite size of the optical fibers of the SDSS multi object 
spectrograph imply that only one quasar in a close pair with sep- 
aration < 55" can be will have a spectrum in the quasar catalog. 
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for the fainter quasar in the pair, appropriate for foUow 
up spectroscopy with a 4m class telescope. 

The left panel of Figure ^ shows the redshift distribu- 
tion of quasars for a flux limit of i < 22.7, i.e. 2.5 mag- 
nitudes fainter than the SDSS flux limit^ of i < 20.2, 
corresponding to magnifications ~ 10. The redshift bins 
used for each source plane are indicated by the shading. 
The right panel of Figure^shows the cumulative number 
magnitude counts of quasars in the redshift bins about 
each of the three source planes. 

4. THE NUMBER OF QUASARS LENSED BY CLUSTERS 

The expected cumulative distribution of lens splittings 
A''(> 9) for the SDSS spectroscopic quasar survey and 
the SDSS photometric survey are shown in Figure[21 The 
green (dotted), blue (dot-dashed), and red (dashed) lines 
show the individual contributions from the source planes 
Zs = 1.0, Zs = 2.0, and Zs = 4.0, respectively. The thick 
solid (black) lines shows the total number of lenses, which 
is just the sum from the three individual source planes. 
The thin solid (black) line shows the total number of 
lenses if the baryons in BCGs are neglected and we ray 
trace through dark matter only. 

BCGs increase the number of lenses by 
^ 50%, consistent with previou s findi ngs of 
dMeneghctti. Bartclmann. fc Moscardini' 120031: ; 

l^l al. Holder, fc Hcnn awi 2004; Ho fc Whit j 12004 : 
IHennawi et al.l l2005a|) . For the SDSS spectroscopic 
quasars, our results are consistent with the existence of 
one wide separation lensed quasar in the sample with 
a splitting 9 ~ 14", namely SDSSJ1004-H4112. This 
lens was discovered in an effective ar ea of ^ 2500 deg^ 
Ijlnada et al.ll2003t lOeuri et al.l |2004D, which is roughly 
half of that used in the left panel of Figure |21 so we 
would predict 1 — 2 such lenses in a quasar sample of 
this size. Our prediction of ^ 4 quasar lenses with 
splittings 9 > 10", suggests that several wide separation 
quasar lenses may still be lurking in the SDSS spectro- 
scopic survey. For the faint photometric quasars, our 
calculation predicts 12 multiply imaged quasars with 
splittings 9 > 10". It is expected that ~ 2 lenses will 
have a separation 9 > 30", and ~ 2 wide separation 
lenses will be at Zg ^ 4. 

At first glance, it seems unexpected that the number of 
quasar lenses in the photometric sample is only a factor 
of 3 higher than the spectroscopic sample, consider- 
ing that the surface density of quasars increases by a 
fact or of ^ 6 and the ar ea is increased by a factor of 
~ 2 ijRichards et al.ll20?)3) . However, recaU that the flux 
limit of the faintest image in the lens, which is identified 
from follow-up observations, is the same for both sam- 
ples, namely i < 21. Most of the lensing probability is 
in the source bins Zg — 2.0 and Zg — 4.0, and the SDSS 
spectroscopic survey flux limit for these bins is already 
i < 20.2, compared to i < 21.0 for Zs = 2.0 and i < 20.5 
for Zs = 4.0 in the photometric sample. Thus the fainter 
flux limits of the photometric sample do not increase the 
number of lenses as much as naively expected. 

In Figure 13 we show the redshift distribution of the 
cluster lenses with 9 > 10" for both surveys. The lens 

^ For the sake of illustration we have used a single flux limit in 
the figure although both quasar samples have different flux limits 
for high and low redshift. 



redshift distribution peaks around = 0.5. Although 
the efficiency for strong lensing peaks at half distance to 
the source plane, which is at z^ = 0.75 for the dominant 
source plane Zg = 2.0, the distribution is skewed toward 
the observer by the cluster mass function. Our results 
for the redshift distribution of the SDSS cluster lenses 
is consistent with the redshift, z = 0.68, of the cluster 
multiply imaging SDSSJ1004+4112. 

5. SUMMARY AND CONCLUSIONS 

We consider the following to be the primary conclu- 
sions of this study: 

1. Our ray tracing calculation through clusters from 
an N-body simulation of the ACDM cosmology pre- 
dicts ~ 4 quasar lenses with splittings 9 > 10" in 
the SDSS spectroscopic quasar sample, consistent 
with the recent discovery of the wide separation 
lens SDSSJ1004-)-4112 which has 9 = 14".6. 

2. Our study predicts ^12 multiply imaged quasars 
with splittings 9 > 10" should b e found in the 
faint photometric quasar sample of IRichards et all 
( 2004). Of these, - 2 should be lenses with 9 > 30" 
and ^ 2 will be at high redshift Zg ~ 4. 

3. BCGs increase the number of multiply im- 
aged quasar lenses by ~ 50%, consistent with 
previous findings in the context of g ia nt arcs 
ilMeneghetti. Bartelm ann fc Mo scardini! I2003bt 
foalal. Holder, & Hcn nawil 12004; .Ho fc Whiti 
l2004HHennawi et al.ll2005a . 

The search for wide separation multiply imaged 
quasar lenses is currently underw ay in both the 
SDSS spectroscopic (Inada et al. 2003t lOguri et alJl2004 
[2005'; Hcnnawi ct al. 2005b ) and photometr ic samples 
lEichiirds et al. 2004; Hcnn awi et alJl2005bD . Looking 
forward, wide field synoptic surveys like Pan-STARRS 
and the LSST will be able to combine variability and 
color selection techniques to detect as many as 10® quasar 
down to i < 27 (Richard Green, private communica- 
tion 2004). Crudely extrapolating from the lensing rate 
computed here, these future surveys are likely to contain 
> 300 quasars multiply imaged by clusters. A compari- 
son of the abundance of these new lenses to theoretical 
predictions as well as detailed modeling of each system, 
will provide a powerful test of the ACDM cosmological 
model. 
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